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A Phenomenological Treatment of Chiral Symmetry Restoration and 
Deconfinement 
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A phenomenological expression for the thermodynamic potential of gluons and quarks is constructed which 
incorporates the features of deconfinement and chiral symmetry restoration known from lattice simulations. The 
thermodynamic potential is a function of the Polyakov loop and chiral condensate expectation values. The gluonic 
sector uses a successful model for pure SU{Nc) gauge theories in which the Polyakov loop eigenvalues are the 
fundamental order parameters for deconfinement. The quark sector is given by a Nambu-Jona-Lasinio model 
in which a constant background Aq field couples the chiral condensate to the Polyakov loop. We consider the 
case of Nf = 2 in detail. For two massless quarks, we find a second order chiral phase transition. Confinement 
effects push the transition to higher temperatures, but the entropy associated with the gluonic sector acts in the 
opposite direction. For light mass quarks, only a rapid crossover occurs. For sufficiently heavy quarks, a first 
order deconfinement transition emerges. This simplest model has one adjustable parameter, which can be set 
from the chiral transition temperature for light quarks. It predicts all thermodynamic quantities as well as the 
behavior of the chiral condensate and the Polyakov loop over a wide range of temperatures. 



Finite temperature QCD, in its variants which 
can be studied with lattice techniques, exhibits 
both deconfinment and chiral symmetry restora- 
tion. Depending on the model under study, there 
may be one or more phase transitions as the tem- 
perature is varied, or no phase transition at all. 
Finite density brings with it an additional rich 
set of possible behaviors at low temperatures. 
The transition from one phase to another may 
be marked by a first- or second-order phase tran- 
sition, or two apparently different phases may be 
connected to one another, as in a liquid-gas tran- 
sition. While we understand much about these 
complex behaviors, we lack a comprehensive, sim- 
ple theoretical tool which gives us insight into 
possible behaviors. In short, we need a Landau- 
Ginzburg theory for finite temperature QCD. 

The behavior of finite temperature QCD can 
be understood from the behavior of two order pa- 
rameters: the Polyakov loop P, associated with 
deconfinement, and the chiral condensate, {ipip}: 
associated with chiral symmetry restoration. Our 
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aim is to construct simple phenomenological mod- 
els incorporating both chiral symmetry restora- 
tion and deconfinement in a unified way, which 
can be applied to both quenched and unquenched 
behavior, for all values of Nc and any number 
of flavors Nf. These models naturally give the 
behavior of P, {'tpip): ^Lud thermodynamic vari- 
able such as pressure as a function of temperature 
and chemical potential. By incorporating both P 
and (V'V') 8-8 order parameters, useful information 
about non-equilibrium behavior is also available. 

Polyakov loop models have been used exten- 
sively to describe the deconfinement transition in 
pure gauge theories ||l| (|] . The two models of the 
deconfinement transition in SU (Nc) gauge theo- 
ries which we have developed 1^. use the eigen- 
values of the Polyakov loop as the fundamental 
order parameters. This accords with the results 
of high-temperature perturbation theory |Q] , and 
our models reproduce the leading T'* behavior of 
perturbation theory at high temperature. Both 
of the models we have studied correctly predict 
a second order phase transition for SU{2) and 
a first order phase transition for SU{Nc) when 
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Nc > 3. 

In the case of SU{3), the center symmetry 
Z(3) aUows us to consider the vacuum expec- 
tation value of the trace of the Polyakov loop 
to lie along the real axis. The trace of the 
Polyakov loop in the fundamental representation 
can be parametrized along the real axis, as 
Pp = 1 + cos(27r/3 — <j>), where <j> — corresponds 
to the confined phase and 4> ^ corresponds to 
a non-zero value for Pp signaling deconfinement. 
We will use our first model, model A, to describe 
deconfinement here. The free energy as a function 
of 6 has the form 
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This model is obtained from a truncation of the 
free energy of massive gluons moving in a con- 
stant Polyakov loop background. Model A has 
a single free parameter, the mass scale M. The 
pressure, energy density and interaction measure 
derived from the above thermodynamic potential 
compare well with pure gauge simulations. 

Quarks break Z{3) symmetry explicitly. This 
effect is small for large quark mass, and the de- 
confinement transition remains first order for suf- 
ficiently large quark masses. Within model A, 
we find that for a single heavy quark, decon- 
finement persists as a first-order phase transition 
until a second-order critical point is reached at 
rric — 2.57M. For quark masses below this value, 
deconfinement occurs as a smooth crossover. 

A complete discussion of QCD at finite tem- 
perature must also incorporate chiral symme- 
try restoration. We include these effects via a 
Nambu-Jona-Lasinio (NJL) model coupled to the 
background Polyakov loop. The free energy of 
Nf = 2 quarks at 1-loop order can be written as 

2 

Vf = ^- NfN,trln{iY'D^ + mo + ct), (2) 

where G is the four-fermion coupling constant, 
the composite field a = and Z?^ is a co- 
variant derivative containing the background 
field. The constituent mass m of the quarks is 
given as m = rriQ + where toq is the current 
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Figure 1. Critical temperature Tc versus M . 
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Figure 2. P and a/lQQMeV versus T in MeV. 

mass. The coupling constant G as well as the 
non-covariant cutoff A required to regularize the 
theory are fixed by zero temperature phenomenol- 
ogy. The values used are [H G = 5.02 GeV~^ and 
A = 653 MeV, which lead to = ^iMeV and 
ilj-ip = —2{2hQMeVY. For sufficiently large tem- 
perature, the finite temperature contribution of 
the functional determinant can be evaluated us- 
ing a high temperature expansion . The quark 
part of our model has no remaining free parame- 
ters once it is calibrated to give the correct zero 
temperature physics. 

The only free parameter in model A coupled to 
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— High-T Expansion 
-- Black Body Limit 
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Figure 3. p/T^ versus T. 

the NJL model is M. The chiral transition tem- 
perature Tc is shown in figure |l| as a function of 
M. As M is lowered, Tc decreases until M = 300 
MeV, below which Tc is essentially frozen at 175 
MeV, the value predicted by the NJL model with 
a trivial Polyakov loop. At very low values of M, 
the energy density divided by changes from a 
monotonically increasing function of temperature 
to one which approaches the blackbody behav- 
ior from above in stark constrast to lattice data. 
If one chooses M — 596 MeV, the value which 
yields a pure gauge deconfinement temperature of 
Td = 270 MeV, then the chiral transition temper- 
ature is 201 MeV which is somewhat higher than 
lattice simulations indicate. Choosing the value 
AI = 300 MeV yields a reasonable value of the 
transition temperature while avoiding the prob- 
lem of poorly behaved thermodynamic functions. 
Once M is fixed there are no free parameters in 
the model. 

The chiral order parameter a and the Polyakov 
loop are plotted as a function of temperature 
in figure Choosing a zero current quark 

mass yields a clear second order chiral transition 
while the deconfinement order parameter displays 
crossover behavior. The pressure divided by 
is plotted in figure ^. The interaction measure 
A = (e — 3p)/T* falls to zero at high tempera- 
tures as l/T^, a behavior also seen in the pure 
gauge theory model. 

The effect of adding a small chemical potential 
can also be studied using a high temperature ex- 
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— High-T Expansion 
-- Extrapolation 
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Figure 4. Tc versus ^. 

pansion of the quark determinant. In figure |4| the 
phase diagram is displayed for M = 300 MeV. 
The chiral transition appears to be second order 
out to at least fi = 130 MeV which is roughly 
the limit of reliability of the high-temperature ex- 
pansion. The dotted line shows an elliptical ex- 
trapolation to higher /i. This line extrapolates to 
precisely the value predicted by the NJL model 
at zero temperature, displayed as a large point in 
the figure. 

We plan to continue studying the phase dia- 
gram of this model. The dependence of the loca- 
tion of the tricritical point on the mass param- 
eter M remains to be worked out. The addi- 
tion of a non-zero bare quark mass is expected 
to change the tricritical point to a second order 
critical point and shift its location slightly in the 
T — fi plane. 
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